Identifying the best PCR enzyme for library amplification in NGS.

Background. PCR amplification is a necessary step in many next-generation sequencing (NGS) library preparation methods [1, 2]. Whilst many PCR enzymes are developed to amplify single targets efficiently, accurately and with specificity, few are developed to meet the challenges imposed by NGS PCR, namely unbiased amplification of a wide range of different sizes and GC content. As a result PCR amplification during NGS library prep often results in bias toward GC neutral and smaller fragments. As NGS has matured, optimized NGS library prep kits and polymerase formulations have emerged and in this study we have tested a wide selection of available enzymes for both short-read Illumina library preparation and long fragment amplification ahead of long-read sequencing.We tested over 20 different hi-fidelity PCR enzymes/NGS amplification mixes on a range of Illumina library templates of varying GC content and composition, and find that both yield and genome coverage uniformity characteristics of the commercially available enzymes varied dramatically. Three enzymes Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier were found to give a consistent performance, over all genomes, that mirrored closely that observed for PCR-free datasets. We also test a range of enzymes for long-read sequencing by amplifying size fractionated S. cerevisiae DNA of average size 21.6 and 13.4 kb, respectively.The enzymes of choice for short-read (Illumina) library fragment amplification are Quantabio RepliQa Hifi Toughmix, Watchmaker Library Amplification Hot Start Master Mix (2X) 'Equinox' and Takara Ex Premier, with RepliQa also being the best performing enzyme from the enzymes tested for long fragment amplification prior to long-read sequencing.

A metagenomic prospective cohort study on gut microbiome composition and clinical infection in small bowel transplantation.

Two-thirds of small-bowel transplantation (SBT) recipients develop bacteremia, with the majority of infections occurring within 3 months post-transplant. Sepsis-related mortality occurs in 31% of patients and is commonly caused by bacteria of gut origin, which are thought to translocate across the implanted organ. Serial post-transplant surveillance endoscopies provide an opportunity to study whether the composition of the ileal and colonic microbiota can predict the emergence as well as the pathogen of subsequent clinical infections in the SBT patient population. Five participants serially underwent aspiration of ileal and colonic bowel effluents at transplantation and during follow-up endoscopy either until death or for up to 3 months post-SBT. We performed whole-metagenome sequencing (WMS) of 40 bowel effluent samples and compared the results with clinical infection episodes. Microbiome composition was concordant between participants and timepoint-matched ileal and colonic samples. Four out of five (4/5) participants had clinically significant infections thought to be of gut origin. Bacterial translocation from the gut was observed in 3/5 patients with bacterial infectious etiologies. In all three cases, the pathogens had demonstrably colonized the gut between 1-10 days prior to invasive clinical infection. Recipients with better outcomes received donor grafts with higher alpha diversity. There was an increase in the number of antimicrobial resistance genes associated with longer hospital stay for all participants. This metagenomic study provides preliminary evidence to support the pathogen translocation hypothesis of gut-origin sepsis in the SBT cohort. Ileal and colonic microbiome compositions were concordant; therefore, fecal metagenomic analysis could be a useful surveillance tool for impeding infection with specific gut-residing pathogens.

NCTC3000: a century of bacterial strain collecting leads to a rich genomic data resource.

The National Collection of Type Cultures (NCTC) was founded on 1 January 1920 in order to fulfil a recognized need for a centralized repository for bacterial and fungal strains within the UK. It is among the longest-established collections of its kind anywhere in the world and today holds approximately 6000 type and reference bacterial strains - many of medical, scientific and veterinary importance - available to academic, health, food and veterinary institutions worldwide. Recently, a collaboration between NCTC, Pacific Biosciences and the Wellcome Sanger Institute established the NCTC3000 project to long-read sequence and assemble the genomes of up to 3000 NCTC strains. Here, at the beginning of the collection's second century, we introduce the resulting NCTC3000 sequence read datasets, genome assemblies and annotations as a unique, historically and scientifically relevant resource for the benefit of the international bacterial research community.

A global resource for genomic predictions of antimicrobial resistance and surveillance of Salmonella Typhi at pathogenwatch.

As whole-genome sequencing capacity becomes increasingly decentralized, there is a growing opportunity for collaboration and the sharing of surveillance data within and between countries to inform typhoid control policies. This vision requires free, community-driven tools that facilitate access to genomic data for public health on a global scale. Here we present the Pathogenwatch scheme for Salmonella enterica serovar Typhi (S. Typhi), a web application enabling the rapid identification of genomic markers of antimicrobial resistance (AMR) and contextualization with public genomic data. We show that the clustering of S. Typhi genomes in Pathogenwatch is comparable to established bioinformatics methods, and that genomic predictions of AMR are highly concordant with phenotypic susceptibility data. We demonstrate the public health utility of Pathogenwatch with examples selected from >4,300 public genomes available in the application. Pathogenwatch provides an intuitive entry point to monitor of the emergence and spread of S. Typhi high risk clones.

Transcriptome of the parasitic flatworm Schistosoma mansoni during intra-mammalian development.

Schistosomes are parasitic blood flukes that survive for many years within the mammalian host vasculature. How the parasites establish a chronic infection in the hostile bloodstream environment, whilst evading the host immune response is poorly understood. The parasite develops morphologically and grows as it migrates to its preferred vascular niche, avoiding or repairing damage from the host immune system. In this study, we investigated temporal changes in gene expression during the intra-mammalian development of Schistosoma mansoni. RNA-seq data were analysed from parasites developing in the lung through to egg-laying mature adult worms, providing a comprehensive picture of in vivo intra-mammalian development. Remarkably, genes involved in signalling pathways, developmental control, and adaptation to oxidative stress were up-regulated in the lung stage. The data also suggested a potential role in immune evasion for a previously uncharacterised gene. This study not only provides a large and comprehensive data resource for the research community, but also reveals new directions for further characterising host-parasite interactions that could ultimately lead to new control strategies for this neglected tropical disease pathogen.

Rapid multi-locus sequence typing direct from uncorrected long reads using Krocus.

Genome sequencing is rapidly being adopted in reference labs and hospitals for bacterial outbreak investigation and diagnostics where time is critical. Seven gene multi-locus sequence typing is a standard tool for broadly classifying samples into sequence types (STs), allowing, in many cases, to rule a sample out of an outbreak, or allowing for general characteristics about a bacterial strain to be inferred. Long-read sequencing technologies, such as from Oxford Nanopore, can produce read data within minutes of an experiment starting, unlike short-read sequencing technologies which require many hours/days. However, the error rates of raw uncorrected long read data are very high. We present Krocus which can predict a ST directly from uncorrected long reads, and which was designed to consume read data as it is produced, providing results in minutes. It is the only tool which can do this from uncorrected long reads. We tested Krocus on over 700 isolates sequenced using long-read sequencing technologies from Pacific Biosciences and Oxford Nanopore. It provides STs for isolates on average within 90 s, with a sensitivity of 94% and specificity of 97% on real sample data, directly from uncorrected raw sequence reads. The software is written in Python and is available under the open source license GNU GPL version 3.

Non-gravitational acceleration in the trajectory of 1I/2017 U1 ('Oumuamua).

'Oumuamua (1I/2017 U1) is the first known object of interstellar origin to have entered the Solar System on an unbound and hyperbolic trajectory with respect to the Sun1. Various physical observations collected during its visit to the Solar System showed that it has an unusually elongated shape and a tumbling rotation state1-4 and that the physical properties of its surface resemble those of cometary nuclei5,6, even though it showed no evidence of cometary activity1,5,7. The motion of all celestial bodies is governed mostly by gravity, but the trajectories of comets can also be affected by non-gravitational forces due to cometary outgassing8. Because non-gravitational accelerations are at least three to four orders of magnitude weaker than gravitational acceleration, the detection of any deviation from a purely gravity-driven trajectory requires high-quality astrometry over a long arc. As a result, non-gravitational effects have been measured on only a limited subset of the small-body population9. Here we report the detection, at 30σ significance, of non-gravitational acceleration in the motion of 'Oumuamua. We analyse imaging data from extensive observations by ground-based and orbiting facilities. This analysis rules out systematic biases and shows that all astrometric data can be described once a non-gravitational component representing a heliocentric radial acceleration proportional to r-2 or r-1 (where r is the heliocentric distance) is included in the model. After ruling out solar-radiation pressure, drag- and friction-like forces, interaction with solar wind for a highly magnetized object, and geometric effects originating from 'Oumuamua potentially being composed of several spatially separated bodies or having a pronounced offset between its photocentre and centre of mass, we find comet-like outgassing to be a physically viable explanation, provided that 'Oumuamua has thermal properties similar to comets.

SeroBA: rapid high-throughput serotyping of Streptococcus pneumoniae from whole genome sequence data.

Streptococcus pneumoniae is responsible for 240 000-460 000 deaths in children under 5 years of age each year. Accurate identification of pneumococcal serotypes is important for tracking the distribution and evolution of serotypes following the introduction of effective vaccines. Recent efforts have been made to infer serotypes directly from genomic data but current software approaches are limited and do not scale well. Here, we introduce a novel method, SeroBA, which uses a k-mer approach. We compare SeroBA against real and simulated data and present results on the concordance and computational performance against a validation dataset, the robustness and scalability when analysing a large dataset, and the impact of varying the depth of coverage on sequence-based serotyping. SeroBA can predict serotypes, by identifying the cps locus, directly from raw whole genome sequencing read data with 98 % concordance using a k-mer-based method, can process 10 000 samples in just over 1 day using a standard server and can call serotypes at a coverage as low as 15-21×. SeroBA is implemented in Python3 and is freely available under an open source GPLv3 licence from: https://github.com/sanger-pathogens/seroba.

PlasmidTron: assembling the cause of phenotypes and genotypes from NGS data.

Increasingly rich metadata are now being linked to samples that have been whole-genome sequenced. However, much of this information is ignored. This is because linking this metadata to genes, or regions of the genome, usually relies on knowing the gene sequence(s) responsible for the particular trait being measured and looking for its presence or absence in that genome. Examples of this would be the spread of antimicrobial resistance genes carried on mobile genetic elements (MGEs). However, although it is possible to routinely identify the resistance gene, identifying the unknown MGE upon which it is carried can be much more difficult if the starting point is short-read whole-genome sequence data. The reason for this is that MGEs are often full of repeats and so assemble poorly, leading to fragmented consensus sequences. Since mobile DNA, which can carry many clinically and ecologically important genes, has a different evolutionary history from the host, its distribution across the host population will, by definition, be independent of the host phylogeny. It is possible to use this phenomenon in a genome-wide association study to identify both the genes associated with the specific trait and also the DNA linked to that gene, for example the flanking sequence of the plasmid vector on which it is encoded, which follows the same patterns of distribution as the marker gene/sequence itself. We present PlasmidTron, which utilizes the phenotypic data normally available in bacterial population studies, such as antibiograms, virulence factors, or geographical information, to identify traits that are likely to be present on DNA that can randomly reassort across defined bacterial populations. It is also possible to use this methodology to associate unknown genes/sequences (e.g. plasmid backbones) with a specific molecular signature or marker (e.g. resistance gene presence or absence) using PlasmidTron. PlasmidTron uses a k-mer-based approach to identify reads associated with a phylogenetically unlinked phenotype. These reads are then assembled de novo to produce contigs in a fast and scalable-to-large manner. PlasmidTron is written in Python 3 and is available under the open source licence GNU GPL3 from https://github.com/sanger-pathogens/plasmidtron.

A brief visit from a red and extremely elongated interstellar asteroid.

None of the approximately 750,000 known asteroids and comets in the Solar System is thought to have originated outside it, despite models of the formation of planetary systems suggesting that orbital migration of giant planets ejects a large fraction of the original planetesimals into interstellar space. The high predicted number density of icy interstellar objects (2.4 × 10-4 per cubic astronomical unit) suggests that some should have been detected, yet hitherto none has been seen. Many decades of asteroid and comet characterization have yielded formation models that explain the mass distribution, chemical abundances and planetary configuration of the Solar System today, but there has been no way of telling whether the Solar System is typical of planetary systems. Here we report observations and analysis of the object 1I/2017 U1 ('Oumuamua) that demonstrate its extrasolar trajectory, and that thus enable comparisons to be made between material from another planetary system and from our own. Our observations during the brief visit by the object to the inner Solar System reveal it to be asteroidal, with no hint of cometary activity despite an approach within 0.25 astronomical units of the Sun. Spectroscopic measurements show that the surface of the object is spectrally red, consistent with comets or organic-rich asteroids that reside within the Solar System. Light-curve observations indicate that the object has an extremely oblong shape, with a length about ten times its width, and a mean radius of about 102 metres assuming an albedo of 0.04. No known objects in the Solar System have such extreme dimensions. The presence of 'Oumuamua in the Solar System suggests that previous estimates of the number density of interstellar objects, based on the assumption that all such objects were cometary, were pessimistically low. Planned upgrades to contemporary asteroid survey instruments and improved data processing techniques are likely to result in the detection of more interstellar objects in the coming years.

ARIBA: rapid antimicrobial resistance genotyping directly from sequencing reads.

Antimicrobial resistance (AMR) is one of the major threats to human and animal health worldwide, yet few high-throughput tools exist to analyse and predict the resistance of a bacterial isolate from sequencing data. Here we present a new tool, ARIBA, that identifies AMR-associated genes and single nucleotide polymorphisms directly from short reads, and generates detailed and customizable output. The accuracy and advantages of ARIBA over other tools are demonstrated on three datasets from Gram-positive and Gram-negative bacteria, with ARIBA outperforming existing methods.

Comparison of classical multi-locus sequence typing software for next-generation sequencing data.

Multi-locus sequence typing (MLST) is a widely used method for categorizing bacteria. Increasingly, MLST is being performed using next-generation sequencing (NGS) data by reference laboratories and for clinical diagnostics. Many software applications have been developed to calculate sequence types from NGS data; however, there has been no comprehensive review to date on these methods. We have compared eight of these applications against real and simulated data, and present results on: (1) the accuracy of each method against traditional typing methods, (2) the performance on real outbreak datasets, (3) the impact of contamination and varying depth of coverage, and (4) the computational resource requirements.

Phylogenetic Analysis of Klebsiella pneumoniae from Hospitalized Children, Pakistan.

Klebsiella pneumoniae shows increasing emergence of multidrug-resistant lineages, including strains resistant to all available antimicrobial drugs. We conducted whole-genome sequencing of 178 highly drug-resistant isolates from a tertiary hospital in Lahore, Pakistan. Phylogenetic analyses to place these isolates into global context demonstrate the expansion of multiple independent lineages, including K. quasipneumoniae.

Inner solar system material discovered in the Oort cloud.

We have observed C/2014 S3 (PANSTARRS), a recently discovered object on a cometary orbit coming from the Oort cloud that is physically similar to an inner main belt rocky S-type asteroid. Recent dynamical models successfully reproduce the key characteristics of our current solar system; some of these models require significant migration of the giant planets, whereas others do not. These models provide different predictions on the presence of rocky material expelled from the inner solar system in the Oort cloud. C/2014 S3 could be the key to verifying these predictions of the migration-based dynamical models. Furthermore, this object displays a very faint, weak level of comet-like activity, five to six orders of magnitude less than that of typical ice-rich comets on similar Orbits coming from the Oort cloud. For the nearly tailless appearance, we are calling C/2014 S3 a Manx object. Various arguments convince us that this activity is produced by sublimation of volatile ice, that is, normal cometary activity. The activity implies that C/2014 S3 has retained a tiny fraction of the water that is expected to be present at its formation distance in the inner solar system. We may be looking at fresh inner solar system Earth-forming material that was ejected from the inner solar system and preserved for billions of years in the Oort cloud.

The TraDIS toolkit: sequencing and analysis for dense transposon mutant libraries.

UNLABELLED: Transposon insertion sequencing is a high-throughput technique for assaying large libraries of otherwise isogenic transposon mutants providing insight into gene essentiality, gene function and genetic interactions. We previously developed the Transposon Directed Insertion Sequencing (TraDIS) protocol for this purpose, which utilizes shearing of genomic DNA followed by specific PCR amplification of transposon-containing fragments and Illumina sequencing. Here we describe an optimized high-yield library preparation and sequencing protocol for TraDIS experiments and a novel software pipeline for analysis of the resulting data. The Bio-Tradis analysis pipeline is implemented as an extensible Perl library which can either be used as is, or as a basis for the development of more advanced analysis tools. This article can serve as a general reference for the application of the TraDIS methodology. AVAILABILITY AND IMPLEMENTATION: The optimized sequencing protocol is included as supplementary information. The Bio-Tradis analysis pipeline is available under a GPL license at https://github.com/sanger-pathogens/Bio-Tradis CONTACT: parkhill@sanger.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments.

Rapidly decreasing genome sequencing costs have led to a proportionate increase in the number of samples used in prokaryotic population studies. Extracting single nucleotide polymorphisms (SNPs) from a large whole genome alignment is now a routine task, but existing tools have failed to scale efficiently with the increased size of studies. These tools are slow, memory inefficient and are installed through non-standard procedures. We present SNP-sites which can rapidly extract SNPs from a multi-FASTA alignment using modest resources and can output results in multiple formats for downstream analysis. SNPs can be extracted from a 8.3 GB alignment file (1842 taxa, 22 618 sites) in 267 seconds using 59 MB of RAM and 1 CPU core, making it feasible to run on modest computers. It is easy to install through the Debian and Homebrew package managers, and has been successfully tested on more than 20 operating systems. SNP-sites is implemented in C and is available under the open source license GNU GPL version 3.

Robust high-throughput prokaryote de novo assembly and improvement pipeline for Illumina data.

The rapidly reducing cost of bacterial genome sequencing has lead to its routine use in large-scale microbial analysis. Though mapping approaches can be used to find differences relative to the reference, many bacteria are subject to constant evolutionary pressures resulting in events such as the loss and gain of mobile genetic elements, horizontal gene transfer through recombination and genomic rearrangements. De novo assembly is the reconstruction of the underlying genome sequence, an essential step to understanding bacterial genome diversity. Here we present a high-throughput bacterial assembly and improvement pipeline that has been used to generate nearly 20 000 annotated draft genome assemblies in public databases. We demonstrate its performance on a public data set of 9404 genomes. We find all the genes used in multi-locus sequence typing schema present in 99.6 % of assembled genomes. When tested on low-, neutral- and high-GC organisms, more than 94 % of genes were present and completely intact. The pipeline has been proven to be scalable and robust with a wide variety of datasets without requiring human intervention. All of the software is available on GitHub under the GNU GPL open source license.

Circlator: automated circularization of genome assemblies using long sequencing reads.

The assembly of DNA sequence data is undergoing a renaissance thanks to emerging technologies capable of producing reads tens of kilobases long. Assembling complete bacterial and small eukaryotic genomes is now possible, but the final step of circularizing sequences remains unsolved. Here we present Circlator, the first tool to automate assembly circularization and produce accurate linear representations of circular sequences. Using Pacific Biosciences and Oxford Nanopore data, Circlator correctly circularized 26 of 27 circularizable sequences, comprising 11 chromosomes and 12 plasmids from bacteria, the apicoplast and mitochondrion of Plasmodium falciparum and a human mitochondrion. Circlator is available at http://sanger-pathogens.github.io/circlator/ .

Roary: rapid large-scale prokaryote pan genome analysis.

UNLABELLED: A typical prokaryote population sequencing study can now consist of hundreds or thousands of isolates. Interrogating these datasets can provide detailed insights into the genetic structure of prokaryotic genomes. We introduce Roary, a tool that rapidly builds large-scale pan genomes, identifying the core and accessory genes. Roary makes construction of the pan genome of thousands of prokaryote samples possible on a standard desktop without compromising on the accuracy of results. Using a single CPU Roary can produce a pan genome consisting of 1000 isolates in 4.5 hours using 13 GB of RAM, with further speedups possible using multiple processors. AVAILABILITY AND IMPLEMENTATION: Roary is implemented in Perl and is freely available under an open source GPLv3 license from http://sanger-pathogens.github.io/Roary CONTACT: roary@sanger.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.